Singlet-Triplet Relaxation in Two-electron Silicon Quantum Dots 
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We investigate the singlet-triplet relaxation process of a two electron silicon quantum dot. In the absence 
of a perpendicular magnetic field, we find that spin-orbit coupling is not the main source of singlet-triplet 
relaxation. Relaxation in this regime occurs mainly via virtual states and is due to nuclear hyperfine coupling. 
In the presence of an external magnetic field perpendicular to the plane of the dot, the spin-orbit coupling is 
important and virtual states are not required. We find that there can be strong anisotropy for different field 
directions: parallel magnetic field can increase substantially the relaxation time due to Zeeman splitting, but 
when the magnetic field is applied perpendicular to the plane, the enhancement of the spin-orbit effect shortens 
the relaxation time. We find the relaxation to be orders of magnitude longer than for GaAs quantum dots, due 
to weaker hyperfine and spin-orbit effects. 

PACS numbers: 72.25.Rb, 03.67.Pp, 68.65.Hb. 85.35.Be 



I. INTRODUCTION 

A promising technology for the implementation of quan- 
tum computation (QC) involves the storage of quantum infor- 
mation in the spin of electrons in quantum dots (QDs). The 
key requirement is that the lifetime of the spins is long com- 
pared with the time required for the operation of logic gates. 
This has motivated the development of dots in silicon, where 
spin-orbit coupling is weak and isotopic enrichment can elim- 
inate hyperfine coupling (HC). Indeed, recent experiments 
demonstrate the presence of long-lived spin states in silicon 
QDs ' . Understanding the processes that relax spins can point 
to strategies for minimizing relaxation and coherence times, 
thereby improving coherent control of quantum systems. In 
the case of electron spins embedded in semiconductor nanos- 
tructures, the relaxation properties are strongly affected by 
the regime of operation. Thus it is important to identify the 
dominant sources of fluctuations in these systems, the mech- 
anisms by which they couple to the spins, and to analyze the 
non-equilibrium decay laws in different regimes of external 
fields. Considerable theoretical work has been performed on 
lifetimes for single-electron spin flip T\ and dephasing T2 for 
GaAs^'^'^ and for Si^. In GaAs these times have been mea- 
sured. Single-spin values for Ti of about 0.5 ms at a field 
of 10 T up to 170 ms at 1.75 T were obtained^^i, while for 
T2 one finds a value of about 1 /js^. A qubit consisting of 
the singlet and triplet states of a two-electron system is also a 
proposal for QC. The singlet-triplet lifetime has been studied 
in GaAs^ '"'''"^^. In particular, Hanson et al. found Ti for the 
singlet-triplet transition in a two-electron GaAs dot to be 2.6 
ms at B=0.02 T. We shall call this Tst- Extensive theoret- 
ical work has been done for Tst in GaAs^ '^'^^''^'^^ and our 
methods are similar to those found in these references. 

In this paper we study the relaxation process for a doubly- 
occupied Si QD in an excited (triplet) state to the ground (sin- 
glet) state, focusing on the computation of Tst. Our main 
motivation is to understand transport through double quantum 
dots. Thus we are concerned with lateral dots defined by gates 
in strained silicon quantum wells. Such dots have a two-fold 
valley degeneracy as well as spin degeneracy, but we will deal 



here with dots where the valley splitting is large compared 
with the first orbital excitation energy. We will focus on natu- 
ral Si with a 4% concentration of ^^Si nuclei, since this is the 
material on which experiments have been performed, but we 
comment on isotopically enriched Si below. 

We assume the levels to be ordered as shown in Fig. [T] The 
relevant energy scales are then the exchange, J = E^i — Ej, and 
the difference between the ground singlet and the first triplet. 
Est = Ej — Es^, where the triplet is formed with a higher en- 
ergy orbital, as depicted in Fig. ^a). The dominant mech- 
anism available in the absence of an external magnetic field 
is the hyperfine coupling with nuclei'**''^ via a virtual stated 
(left arrows of Fig. [Db). HC cannot cause a direct T ^ S 
transition because the nuclei cannot absorb the energy. So 
the change in energy of the electron spin must be accompanied 
by the emission of a phonon^Si^ii^. The electron-phonon in- 
teraction mixes thus different orbital states via a deformation 
potential in this process, while the spin-flip is provided by the 
HC. This is the dominant process at zero applied magnetic 
field. A second relaxation channel is through spin-orbit (SO) 
coupling. SO coupling mixes different spin states through the 
Rashba SO coupling^^. This leads to a non-vanishing matrix 
element for the phonon-assisted transition between a singlet 
and a triplet state in the absence of time- reversal symmetry 
(right of Fig. [T]! leading to an increase in the relaxation rate. 
Est as the field is increased. Our aim here is to compute the 
singlet-triplet relaxation rates due to these two mechanisms as 
a function of external field. 

We outline the method in the following section, and justify 
the approximations that are made. We then present results 
and discussion. 



II. METHOD AND RESULTS 

We consider first the case of low fields. For reasons to 
be discussed below, we may neglect the spin-orbit coupling 
and the Hamiltonian is written as T/q + 5//, where Hq con- 
tains the confining potential of the QD and hH = T/hc +^^ph- 
Here T/ph corresponds to the electron-phonon coupling (which 
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FIG. I', (a) Level scheme for a doubly-occupied silicon QD: the ground sin- 
glet involves only one orbital, whereas an additional orbital is needed to form 
the triplet and higher states, increasing the level separation, (b) Dominant 
processes in the quantum dot. The two rates indicated are the combined 
hyperfne and phonon rate, r^c-ph <^nd the combined spin-orbit and phonon 
rate, Ty^^-^h- The energy separation of the ground-state singlet and first 
triplet is denoted by Est = -Et — ^^s., and the exchange splitting by J. 



conserves spin), and i/hc is the hyperfine coupling that causes 
spin mixing in the dot. The confining potential is taken as 
parabolic with circular symmetry. We do not consider the 
Coulomb interaction explicitly in the Hamiltonian, as was 
done, for example, by Golovach et al.^, whose main interest 
was in the regime close to the singlet-triplet crossing. Instead 
the interaction is included phenomenologically through a pa- 
rameter J, the singlet-triplet splitting. This is a crude form of 
mean field theory, but is reasonable as long as we are far from 
the singlet-triplet crossing point. Far from this point there are 
no energy denominators that depend sensitively on the inter- 
action strength, and matrix elements depend smoothly on the 
strength. 

For the purpose of this paper (QD formed in a bi-axially 
strained quantum well grown along the z-axis), we consider 
the lowest electric subband, so that the wavefunction in 
the z-direction is fixed. We utilize the Fock-Darwin (FD) 
states for the lateral dimensions (j)'" ™) ix,y) to construct our 
wavefunctions, {'P,}, that diagonalize Hq. For these states, 
orbital and spin degrees of freedom factorize: \^Sg) = X{z) x 
[|1)|1>] ® \S), \^>t) = Xiz) X [|1>|2) - m)]/V2® \T), and 
\W)=x{z)x[\l)\2) + m)]/V2^\S). Here, |5) ^ [| U 
)- iT)]/V2 and \T+fi-) ^ [\ TT);(| U)+ iT»/V2;| U); 
denote spin states, |1) = ^^^'^''{r) and |2) = (|)('^'=^')(r), so 
|1)|2) = (|)(°'°)(ri)(|)(''^±')(r2). The electron-electi'on inter- 
action is taken into account only phenomenologically through 
the parameter J. The amplitude of a transition between the 
triplet 1^7) and singlet \^Sg) ground state via an excited 
state is given in second order perturbation theory. 
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The transition rate from \^t) to \^Sg) is then given by Fermi's 
golden rule: 
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In this notation, denotes the initial state of electron, nu- 
clei andphonons, = l^j-) (g) \i„) ® |/ph), likewise = 
\^Sg) ® l/n) l/ph)- The av subscript indicates that the initial 



states of the nuclear and phonon systems are averaged over 
thermal ensembles, and that the final states of these systems 
are summed over. In this paper we take the temperature to 
be 100 mK, as this is roughly the temperature at which ex- 
periments are done. The chief approximations involved in 
the calculation are the use of second-order perturbation the- 
ory and the truncation of the Hilbert space to just two singlet 
states and one triplet state. The first approximation is ex- 
cellent - the rates turn out to be on the order of seconds; at 
those time scales the Born-Markov approximation implicit in 
Golden-Rule calculations is surely valid - the time scales in 
the bath are probably of the order of the time for a phonon to 
traverse the dot. The validity of the second approximation is 
less clear - in high-symmetry dots such as we are considering 
here the phonons do not couple to highly excited states in the 
dipole approximation, but real dots may be more disordered. 

For silicon under compressive stress along [001], the elec- 
tron interacts with a phonon of momentum q via deformation 
potentials^-^"*-^^ so the Hamiltonian reads: 
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{Zj + Zu)e%), (3) 



where (n^y — 1 jfl^^ln^) = {hng/lMc&q), Mc is the mass of the 
unit cell, Uq is the phonon occupation number and h(£)q is the 
phonon energy. Here s denotes the polarisation of the phonon 
(two transverse and one longitudinal), q is the wavevector, 
and E„ and are the electron-phonon coupling parameters. 
This is slightly simpler than the corresponding Hamiltonian in 
GaAs because of the absence of the piezoelectric coupUng in 
(centrosymmetric) Si. 

Next, we evaluate the spin-flip matrix element given by 
(^s/|//hc|1'7'), which is provided, at low magnetic fields, by 
the bath of nuclear spins of the "^Si isotope'**. Accordingly, 
we consider a contact Hamiltonian, 

H^o=Z'^liBmS,lMr,-Rj) =A'£SiIj5{n-Rj), (4) 
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where 5, (Ij) and r; (/?,)denote the spin and position of the ith 
electron (7th nuclei), and r| and A are hyperfine coupling con- 
stants. Inserting eq. (|4|l and (|3]l into (|2|l, we get an expression 
for the singlet- triplet rate: 



TsT = Fph X 



Y,[\{T\[srir5{r,-R,)]\s')\\ (5) 



Fph describes the phonon rate derived from (|3]l, and mixes the 
different orbitals contained in ji') and \Sg). We use the electric 
dipole approximation (e*'*'' w 1 ± iqr) in[3] which is valid for 
the range of energies considered here (~ 200/jeV— ), 
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2ps,(2Jt)2'^i Jo (£>q 
+|(/|z|5,)|2(e-)2]5(nco-(7 + esT)). (6) 
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To evaluate the integral over momenta of (O, we assume an 
isotropic phonon spectrum, Ep^ ~ ftoj^s and a linear dispersion 
relation, cOfy y = Vj^', being the sound velocity of the mode s. 

The sum over j of ^ can be transformed to an integral by 
introducing C„ as the volume density of ^^Si nuclei, resulting 
in a compact expression for triplet-singlet relaxation: 



TsT 



(7) 



where y^., contains the result of the angular integral 
which depends on the mode s along the coordinate 
i- lix = Jiy = 47t(E2/3 + 2E,,E„/15 + ^2/35), y,, = 
4% (s2/3 + ZA/5 + Zl/7), y,^ = y,,, = 4kZI/35, and y,^. = 
47tS2/70. 

It is important to note that Est is proportional to C„, i.e., 
to the total number of nuclei N„ with which the electrons in- 
teract. This is consistent with the simple picture that the re- 
laxation rate is proportional to the mean-square fluctuations in 
the random hyperfine field. Formulas for spin relaxation rates 
due to hyperfine coupling that give an apparent proportional- 

-1/2 

ity to N„ are common in the literature, and have given rise 
to the incorrect notion that some sort of motional narrowing 
is at work. This is not possible, since the fluctuations in the 
nuclear spin system are slow. In any case the rate must van- 
ish as A^„ — > 0. These formulas are correct, but they generally 
involve other parameters that actually vary with N„ . 

Fig. |2]represents Tst — obtained as described in Eq. |7] 
as a function of the level separation, J = E^i — Ej for a given 
Est- In case of Si, r] = 186,^^ yielding A « 2 x lO^^eV-nm^. 
Only about 4% of the nuclei have spin, so C„ « 0.04 x 8 /vq ~ 
2 nm^^ (v'o « 0.17 nm^). Other parameters used are: E„ = 
9.29 eV, = -10.7 eV, p,. = 2330 kgm-^ v, = 9330 ms"', 
Vf =5420 ms^^ Est = 200 /jeV. At small J, Tst increases 
as a function of J : the triplet |7± o)and singlet \S') levels are 
strongly mixed by the hyperfine interaction and the phonon 
density of states increases as a function of level separations. In 
the limit 7=0, the rate is given by phonon emission, which we 
found to be of the order of lO'^ s^^ Thus Tst appears very 
small as y ^ in Fig. |2] but it is not zero. At large J, spin 
mixing is lessened because there is an energy denominator and 
phonon relaxation is then suppressed by spin conservation. 

The calculations are obviously consistent with the observed 
lower bound of Tst > 15 fjs given in Ref. ^ . It is expected that 
more stringent bounds, hopefully also upper bounds, will be 
available soon. 

We now move to the case of finite applied magnetic field 
B. We first take the the field along the growth direction 
(perpendicular to the 2-DEG): B — B±. This allows a di- 
rect r — > 5 transition to occur in the presence of a Rashba 
SO coupling^**. The Rashba field is a consequence of struc- 
tural inversion asymmetry^^ in the heterostructure. Note that 
no bulk inversion asymmetry needs to be considered in a cen- 
trosymmetric crystal like Si. The SO Hamiltonian due to the 
full confining potential for the device considered here where 
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FIG. 2: Tst a function of J for the process depicted in the inset: when 
levels \T) and \s') are close in energy (J — > 0), the process is dominated by 
hyperfine coupling (fast), and as this energy splits, the process is dominated 
by the phonon emission, . For higher energies, the dipole approximation 
breaks dowr& 



la/rQD ^ 1, 'a being the Data-Das device length^", is then 
given by^ 

- a(a X p)-. 



(8) 



H mixes the spin states \Sg) and \T ), and the orbital wave- 
functions as well, so virtual transition to a higher state \s') is 
no longer needed and we have: 
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^^ = ±g^iBB. 



EST + A± 

(9) 

Here, the phonon just ensures energy conservation. Notice 
that (^glS/Zsol?^) is zero if B = 0. This 'Van Vleck cancella- 
tion' has been known for decades^i*^, but has been clarified 
in recent years^^ particularly by performing a spin-dependent 
unitary transformation in which the first-order term in a is 
eliminated^^. In explicit calculations in the original basis the 
cancellation occurs due to the fact that the admixture of r+ 
and T is equivalent in magnitude but opposite in sign. The 
key point is that spin-orbit-induced transition rates are always 
proportional to B^ (or higher powers of B in the case of spin 
l/2dots)^ 

Proceeding as in Eqs. (1) to (|5]i, we find the rate for the 
direct transition: 
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Note that for B = F^j vanishes. Since the phonon does not 
mix different orbitals in zeroth order in the multipole expan- 
sion to zeroth order (i.e., e'^'^'' « 1), Fpj^ reads now: 
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{Z,, + Z„)eiq,)f8{fm-esT) (10) 

Inserting (|9]l into dU and evaluating ( fTOl i. we get the spin- 
orbit-induced relaxation rate, which is quadratic in B and a, 
in agreement with earlier treatments^: 
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with Y, = 47I(S2 + + Zl/5) and = 47t(E2/5 + 

4S,,E„/15 + 22^/15). Here we use a « 50 m/s, fol- 
lowing Ref.^ and adapting the result to an electric field of 
lO^^V/m, an estimate for for a QD with a 2DEG density 
of ~4x 10^'cm^^. (One should note, however, that this value 
of a is very uncertain.) Fig 3 (a) contains in black the relax- 
ation in the absence of the Rashba coupling (a = 0) and in red 
the relaxation with the additional Fg-j. (Fst + Fg.j.) . We plot 
the results for two values of J. For a = 0, Ist is a relatively 
weak function of 7, since the hyperfine interaction is not very 
sensitive to field. At finite a, Bx activates the mixing and Tst 
decreases rapidly. 

If the magnetic field B|| is parallel to the 2DEG, the only 
effect is that the spin splitting increases, and larger relaxation 
times are obtained (Fig. ^p), as this increases, on average, 
the energy separation E^j±^ — E^^') (see diagram). In contrast, 
the perpendicular magnetic field Z?x decreases the relaxation 
time, as we have seen. This anisotropy in applied field of 
Tst would be a critical signature of the spin-orbit effect. We 
compare our result for a B|| = 0.02 T (Tst ~ 500 ms) to the 
experimental value obtained for GaAs (2.6 ms, '") and find that 
Si has a singlet-triplet relaxation time more than two orders of 
magnitude larger than GaAs. 
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FIG. 3: (Color online) Tst as a function of B for two values of J: The solid 
lines correspond to a plausible value ofJ= 0.02 meY and the broken lines to J 
= Est =0.2 meVfor comparision, (a) as a function ofB\\ and (b) as a function 
ofBj^, including Rashba (red) and without Rashba (black). Tst increases 
with B|| and decreases with B^. 

Note that the behaviour of Tst here is different to the 
much studied GaAs-based devices^ii^ii^ because of the nature 
of spin-orbit coupling and electron-phonon coupling in non- 
centrosymmetric materials: The BIA is absent and there are 
no piezo-phonons, also avoided crossings of the singlet and 
triplet energy level does not occur for the magnetic fields con- 
sidered here, giving a monotonous behaviour. 



III. DISCUSSION 



There are several ways to measure Fst. In single-dot sys- 
tems, this can be realized using a single pulse^. Alternatively, 
one may use the following sequence: in the first phase, the 
state can be prepared so that only one electron is present in the 
QD, and in the next phase, the triplet would be available for 
conductance, unless it relaxes to the singlet. Measurement of 
the current for different values of the pulse duration then gives 
a direct method to determine Fst. The latter experiment has 
been performed in GaAs^^. In double-dot experiments, Fst 
is one of the parameters in the rate equations that determine 
the measured current, so these experiments also provide an 
avenue for the determination of the singlet-triplet lifetimoii. 

One should note immediately that Tst is considerably 
longer in natural Si than in GaAs, generally by orders of mag- 
nitude. This is expected in a system with weaker spin-orbit 
coupling and fewer spinful nuclei. The times we find are of 
the order of seconds for the most part. It is possible to reduce 
the time by applying a perpendicular field, which can serve as 
a very useful diagnostic. It is also possible to lengthen Tst by 
the use of isotopically enriched Si, i.e., pure ^*^Si. This would 
eliminate the hyperfine mechanism but it would not get rid 
of spin relaxation entirely, as higher-order effects of SOC are 
still present even at B = 0. However, these effects are quite 
small in Si. It seems likely that other effects such as flux noise 
will be the limiting factor in isotopically enriched Si. 

In summary, we have calculated the dominant rates for 
phonon-assisted triplet-singlet relaxation of a silicon quan- 
tum dot. Tst is found to be of the order of hundreds of ms, 
very sensitive to the exchange energy 7, and even longer in the 
presence of a B||, to seconds. In the presence of a Bx> ^ di- 
rect transition becomes possible, increasing (decreasing) Fst 
(Tst). Due to weak spin-orbit and hyperfine coupling, silicon 
offers very long coherence times, which are required for solid 
state qubits. 
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